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Quantum Weak Energy Inequalities (QWEIs) have been established for a variety of quantum 
field theories in both flat and curved spacetimes. Dirac flelds are known (by a result of Fewster and 
Verch) to satisfy QWEIs under very general circumstances. However this result does not provide 
an explicit formula for the QWEI bound, so its magnitude has not previously been determined. 
In this paper we present a new and explicit QWEI bound for Dirac fields of mass m ^ in four- 
dimensional Minkowski space. We follow the methods employed by Fewster and Eveson for the 
scalar field, modified to take account of anticommutation relations. A key ingredient is an identity 
for Fourier transforms established by Fewster and Verch. We also compare our QWEI with those 
previously obtained for scalar and spin-1 fields. 

^ . PACS numbers: 03.70.+k, ll.lO.Ef 
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D ■ I. INTRODUCTION 
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' It has been known for a long time that the energy density of a quantum field may assume negative values 1] . If the 
' . magnitude and the duration of such negative energy densities were unconstrained various exotic phenomema might 
■ be possible, ranging from the manufacture of macroscopic traversable wormholes to large-scale violations of the 
^ , second law of thermodynamics Q . 

QQ ■ However, quantum field theory itself contains mechanisms to limit the extent of negative energy densities: certain 
^\ ' weighted averages of the stress-energy tensor obey state-independent lower bounds known as quantum weak energy 
' inequalities (QWEI) [or simply quantum inequalities (QI)] which severely limit exotic phenomena such as those 
C ■ described above. 

. Bounds of this type have been obtained by various means 0, IE IE S H IS Qi Oi Q for the scalar, 

electromagnetic and Proca fields, in both flat and curved spacetimes, leading to results of great generality. Taking 
four-dimensional Minkowski space as a specific example, and averaging the energy density along an inertial worldline. 
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^ the QWEI satisfied by these three theories may be written in the form 8, 9, 11, 12] 

J dt{:no:)4t,xo)g{tf J du\g{u)\^u*Q:,{u/m) (1.1) 

for any smooth, real- valued g vanishing outside a compact region, where m is the particle's mass, S denotes its helicity 
(5 = 1 for scalars, 2 for photons and 3 for massive spin-1 particles), while 

S ■ 1/2 

Q3(x)=fl-^') (l-^)-J-\n{x + V^^) (1.2) 



and is replaced by lima;_»oo Qsix) = 1 in the massless case. (In curved spacetimes the QWEIs of these theories are 
not simply related by an overall factor.) 

In the spin-i case only two arguments are known. One, due to VoUick (based on earlier work of Flanagan |E]), 
exploits the conformal properties of massless Dirac fields in two-dimensional spacetimes to obtain a sharp (and explicit) 
QWEI for an arbitrary sampling function. The other, due to Fewster and Verch [l6|, establishes the existence of 
QEWIs for Dirac and Majorana fields of arbitrary mass in general four-dimensional globally hyperbolic spacetimes by 
drawing on techniques from microlocal analysis. However, this approach did not lead to an explicit formula for the 
QWEI bound, so its magnitude — even in four-dimensional Minkowski space — has not previously been determined [20j |. 
(Specific classes of Dirac states with local negative energy densities have been studied 0, and shown, in the 
massless limit, to satisfy the scalar field QWEI obtained in of course, this does not amount to a proof of the 
QWEIs.) Thus, for example, it has not previously been clear whether or not it is more or less easy to support a 
traversable wormhole using Dirac fields, rather than with the scalar or electromagnetic fields. 
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In this paper, we present a new QWEI derivation for Dirac fields, which combines ideas taken from with 
an approach first developed for the scalar field in 0, 0. The argument is straightforward, involving the use of 
anticommutation relations and an identity for Fourier transforms, and yields an explicit, closed-form, bound. In 
four-dimensional Minkowski space this bound is 

1 

dt{:Too:)^{t,Xo)gitf ^^—^ / du \g{u)\^ Qi{u/m), (1.3) 



where 



3/2 

(x) = 4 I 1 - - ) - 3Q3(x) (1.4) 



is replaced by limx^oo ^3(2:) = 1 in the massless case. Although we have restricted ourselves, for simplicity, to 
Minkowski space, our technique certainly extends to static spacetimes and may even suggest a new approach in the 
general globally hyperbolic case. It is hoped to return to these issues elsewhere. 

Conventions: We take the metric to have signature (H ) and adopt units in which ft = c = 1. The Fourier 

transform / of a function / on R is defined by 

/•oo 

/(^)= / dt/(t)e— *. (1.5) 



II. THE QWEI DERIVATION 

To establish our notation, we begin by recalling some basic facts about the Dirac fields. The Dirac equation for a 
fermion field ip of mass m is 

(^7*^5^-771)^ = 0, (2.1) 
where the 7-matrices are given in terms of the Pauli matrices (7^ by 

7»^(j_°i). „■) (.^1,2,3) (2^2) 

and obey {7^,7"^} — 277'"'. The stress-energy tensor is 

i < — > < — 5- 

T^,^ = [V'7m dvi^^ i^lv df, tp] , (2.3) 

where 7/) = so the classical Dirac energy density is 

no ^ '-[i^^ - i'H]. (2.4) 

(Of course, this quantity is unbounded both from above and below, in contrast to the energy density of a classical 
scalar field. This is at the root of the differences between the QWEIs obeyed by the scalar and Dirac fields.) To 
simplify our notation, we quantise the Dirac field in a box of side L and then take the continuum limit L —^ 00 towards 
the end of the argument. We emphasise that this is done purely for notational convenience. The Dirac field operator 
ipix) is therefore defined in terms of creation and annihilation operators by 

^(^) = E E [fea(fc)7."(fc)e-'^-+4(fc)7;"(fe)e*-], (2.5) 

k Q = l,2 



where {L/Ti)k runs through 1? , and the four-vector k has components fc" = {ujk,k) with ujk — -t- tti^. The 

annihilation and creation operators satisfy the anticommutation relations 

{6„(fc),6],,(fe')} = <5„,a'<5fc,fc-l (2.6) 
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and 



(2.7) 



with all other anticommutators vanishing. Finally, the spinors u°'{k) and where a = 1,2 labels the two 

independent spin states, are given by 



/ 



ii"(fe) 



cr ■ k 



(. 



CT ■ k 



and w"(fe) 



\i/2wfc(c^fc + ra)V j 



V 



I 



(2.8) 



in which the two dimensional column vectors 0" are (/>^^ = (1, 0) and = (0, 1) and V — . The normalisation has 
been chosen so that 



^|lu"(fc)||2=^||«"(fe)|p^-. 



(2.9) 



Substituting the quantum field ip(x) into the classical energy density Eq. H2.4(l and normal ordering according to 
the prescription 

:d„(fc)4(fc'): = -4'Wa(fc), 



(2.10) 



we obtain the normal ordered energy density at the spatial origin (t, 0) in the form 

:roo:(i,0) = ^ 5] 5]{(c.,+c.;,0[&l(fe)^^a'(fe')^^"nfcK'(fc')e'("'= 

k,k' a^a' 

+ dl, {k')d^{k)v^^ {k)v"' (fe')e-^("'»-"^')*] 

+ (cc-fc, - Wfc)[d„(fc)&„,(fc')«"^(fc)«"'(fe')e"'^"'=-^"^'^* 
-6t(fe)4.(fc')""^(fc)«"'(fc')e'^"'=+"'='^*]}- 
We will consider weighted energy density averages of the form 



Tf = 



dt:Too:it,0)fit), 



(2.11) 
(2.12) 



measured by a stationary observer at spacetime origin 0, where / is a non-negative sampling function. Thus, 

k,k' a.a' 

+ dl{k')do.{k)v"\k)v"' {k')f(iUk - UJ^')] 

+ {lO^, - C^fe) K {k)ba' {k')v"^ {k)u''' {k')f{LOk + t^fc') 

~bi{k)dl,{k')u'^\k)v-\k')j{-ujk-u:^')]}. (2.13) 

Our aim is to determine a lower bound on expectation values of Tf in the case where / = .9^, for some real- valued, 
smooth, compactly supported function g. To this end, we first define a family {O^ij/i e R, i = 1, 2, 3, 4} of Fock space 
operators by 



= J2{d{-L0k' + ^^)bo.'{k')uf (k') +giu;^, + ^l)dl{k')vf (k')}. 

a' 

A short calculation, using the anticommutation relations and the normalisation H2.9|) . shows that 
OlO^, = S^l + E E {Si-^'^ + M)5(-^fc' + M)6l(fc)&a'(fc')«"t(fc)«">') 

- g{LUk + fi)g{ujk' + fi)dl, {k')d^ (fe)«"t (fc)^;"' (fc') 
+ g{uk + y)g{-uJk' + y)d^ {k)b^' (fc')«"^ (fc)^"' (fc') 
+g{-u:k + ^AK^k' + P)bi{k)dl,{k')u'^\k)v'^' {k')] 



(2.14) 



a,a' k.k' 



(2.15) 
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where 

V 



S^. = ^Y.\9(^>^+^')\^^ (2-16) 



k 



and we have imphcitly summed over the spinor index i. 

To relate this to our expression for TV we use the following Lemma, which plays an analogous role here to that 
played by the convolution theorem in 0, Q • 

Lemma 1: If f ^ for some real-valued, smooth, compactly supported g G C^(M) then 



{uj + uj')f{uj - u') ^ — Hg{uj - fj.)g{uj' - fj,). (2.17) 



This Lemma was also used in the QWEI derivation in [T^: for completeness we give a proof of this statement in 
the Appendix which differs slightly from that given in that reference. Using this Lemma with the anticommutation 
relations Eq. (|2.6|l and Eq. H2.7() . we calculate 

/OG ^ ^ 

dfifi{OlO^-S^l) = 7r^^{K+c.;,0[&l(fc)V(fc'K^(fcK'(fc')7(^fc'-^fc) 

+ 4,(fc')da(fc)""^(fc)«"'(fc')/K - ^k')] 

+ {LJ,,,^UJk)[dJk)b^'{k')v''Hk)u''\k')f{LUk+LU^') 

-bUk)dl{k')u-^ik)v"'ik')fi^u;k-u;^,)]}. (2.18) 

Comparing with Eq. H2.13|l , we conclude that 

Tf = — d^i^i{olo^~s^l), (2.19) 

which is similar in form to the decomposition of the scalar field energy density employed in @ . Next we compute the 
anticommutator 

{ol,o,.} = Yl T.iah^'^ + i^)d{-^k' + ii)mk),b^ik')}u^\k)u^' {k') 

a,a' k.k' 

+ giiUk + Ai)?Kr:M{rfa(fe), 4, {k')}v"\k)v''' (fc')} (2.20) 

where we have again implicitly summed over the spinor index i and used the fact that anticommutators between the 
&a(fc) and da(k) vanish. Using the anticommutation relations, the above expression reduces to 

{Ol,0,,} = j^|?(-c.,+^)|2|l^"(fe)|l2 + ^|5(c.;, + /i)nl^;"(fe)|M 1 

= (^-p + 5^)l, (2.21) 

where we have also used Eq. H2.9|l and the fact that \g{u)\ is even (because g is real). Splitting the /i-integral of 
Eq. H2.19|l into two pieces and using Eq. (|2.21|) gives 

^f = ^J^ i^ioUo^. - 5^1) + — d^i^i[{s^ + 5_^)i - o^.,ol - s^i] 

1 1 

= 2^i, diJiiJi{olp^,-s^t) + — d^i^lis^^l-o,,ol). (2.22) 

Now, for any yit > 0, we have fi{Oj^^0^i)jp ^ 0, while for any fi ^ 0, we have — ^{O ^iO'^^j) ^ ^ for all physically 
reasonable states |23|- Thus, 

I If" 
i'^f),p > dl^l^S^ + — J dfifiS^^, 



1 



dnnS^, (2.23) 
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so we have obtained a state-independent lower bound on {Tj)^. We now proceed to calculate the right-hand side of 
Eq. 12.23|l and show that it is finite. Inserting the definition Eq. H2.16|l of 5*^ into Eq. (|2.23|) . we find 



2 1 

TT Jo V ^ 



(2.24) 



Taking the continuum limit 



the inequality becomes 



2 noc nCX} 



(27r)^ 

3 />OC 

d/i/i / dijjuj\/ iJ^ — w?-\g{uj -I- /i)p, (2.26) 



2 

TT 

_ 1 

where we have made a change of variables in the integration to spherical polar coordinates, integrated over the angular 
variables and changed variable again to the energy — \/k?' + mn?. The next step is to make a further change of 
variable 

u = uj + ^i, V = Lo, (2.27) 
after which the quantum inequality Eq. H2.25|l becomes 

(Tf) ^ I du\g{u)\'^ I dvv{u — v)\/v'^ — m? 

(2.28) 



4,/;<,„BWp(i.,(.,'-™V'^-i.«<53(^)), 



where the function Qz{x) is given in Eq. p.2|l . 

Using the translational invariance of the theory to move from a worldline (t, 0) to (t, a;o), we have therefore 
established the bound 

dt{:Too:)4t,Xo)g{tf^-j^J^ du\g{u)\^u*Q!^ {u/m), (2.29) 

for any Xq^ where 

Of (x) = 4 fl - 4 V - 3g3(x), (2.30) 



as claimed in the Introduction. This bound is clearly finite for any smooth compactly supported g, as u^Q^ {u/m) 
grows like for large u, while |5(u)p decays faster than any inverse polynomial in u as m ^ cxd. 



III. MASSLESS CASE 



We now look at the simplest case of the general QWEI derived above; the massless case. We then compare our 
result for the Dirac field with those for the scalar field derived by Ford and Roman 0, |^ and Fewster and Eveson |^ . 
Substituting to = in the first line of Eq. H2.28() and performing the integral over v, we immediately find that the 
four-dimensional QWEI for massless fields takes the form 

1 f°° 

(Tf)^ > du\g{u)\\\ (3.1) 
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FIG. 1: Plot of 



which is also the m ^ hmit of Eq. H2.29|l because Q!^{x) — > 1 as x ^ cx). Since the integrand is an even function in 
u, we may extend the range of integration to the whole of M to give 



1 



247r3 
1 

12^ 



— oo 
oo 



dt\g"{t)\' 



(3.2) 



where we have used Parseval's theorem to yield a t-space version of the quantum inequality in the last line. In passing, 
we observe that, because is bounded between and 1 (see Fig.nj, the introduction of a mass only serves to tighten 
the QWEI bound. Accordingly, Eq. (|3.2|l also holds for massive Dirac fields. 

Comparison with the quantum inequality for the scalar and electromagnetic fields in Minkowski space derived 
0, ^3 shows that the massless Dirac bound is weaker than the corresponding scalar bound by a factor of | , and 



m 



4 
3' 

stronger than the electromagnetic bound by a factor of |. Since none of these bounds are optimal, it is difficult to 
draw definite conclusions from this beyond the observation that the bounds are of comparable magnitude. 
In particular, for the choice of a Lorentzian sampling function peaked at < = 0, with characteristic width r 



in which t sets the averaging timescale, it is straightforward to calculate 

+2 ^2\2 



1 



127r2 



dt 



{2t' - t'Y 
7r(t2 + r2)^ 



36 



20487r2r4 



(3.3) 



(3.4) 



(at least for those states for which the right-hand side is well-defined). This is a slightly stronger result, by a factor of 
than the quantum inequality derived by Ford and Roman 0,0] for the scalar field. We have therefore explained why 
the specific classes of states considered in ^3,^3 satisfy the Ford-Roman QWEI. We also observe that a traversable 
wormhole supported by a Dirac field would be subject to slightly tighter constraints than those envisaged in 0. 
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IV. CONCLUSION 



We have derived a new and explicit QWEI for the Dirac field of mass m ^ in four-dimensional Minkowski space, 
by straightforward means. Our argument is related to the general (but inexplicit) QWEI obtained in in that 
Lemma 1 is at its heart; one may therefore regard it as an improved and streamlined version of |16| (for the particular 
case of averaging along an inertial trajectory in Minkowski space) with the added benefit of a closed-form expression 
for the bound. Although it is at present somewhat formal, we expect that our argument can be made fully rigorous 
(by means of microfocal techniques, cf. |lft IT^ ) and that it can shed further light on the general situation in curved 
spacetimes, again leading to explicit bounds. 

Our argument is also related to that developed for the scalar field in and which we briefiy summarise 

in Appendix IbI In each case, an identity for Fourier transforms (the convolution theorem or Lemma 1) and the 
(anti)commutation relations are used to express the averaged energy density as an integral of manifestly positive 
operators, modulo a c-number term which eventually provides the lower bound. Extra care is required in the Dirac 
case because the classical energy density is unbounded both from above and below. 

Finally, in contrast to Vollick's two-dimensional bound .15], our QWEI is not expected to be the sharpest bound 
possible. Again, it is hoped to address this issue elsewhere. 



APPENDIX A: PROOF OF LEMMA 1 



In this appendix, we wish to state and prove the following lemma (originally proved in by a slightly different 
method): If / = for real- valued g € C^(M), then 



(A + A')/(A - A') = / -^figiX ~ m)5(A' - m)- (Al) 



To obtain this result, first notice that the left-hand side is the momentum space integral kernel for the operator 

A=pf + fp (A2) 
where / acts by multiplication and p = —id/dx. That is, 

1"°° d\' ^ 

J-oo 

holds (at least) for all smooth, compactly supported ip. On the other hand, because f — g^, we also have 

A = pg^ + g'^p = 2gpg + [p, g]g + g[g,p] = 2gpg (A4) 
because [g,p] — ig' , another multiplication operator. Written as an integral operator, we therefore have 

/?(A-m)m/ —g[^Ji~X'm\'). (A5) 

Since both ljA3|l and IjASp hold for all smooth compactly supported -0, the two integral kernels must also agree, 
therefore yielding \K\\ . using the fact that g(^ — A') — g{X' — ji) because g is real. 



APPENDIX B: SKETCH OF THE SCALAR FIELD ARGUMENT 



For ease of comparison with our Dirac QWEI, we briefly summarise the scalar field argument of |8|,|9j. (Notation 
differs slightly, and we have specialised to four-dimensional Minkowski space.) This argument begins by observing 
that the averaged energy density may be written as a sum of terms of the form 

± f{uok + ioMk')p{k)a{k')a{k) + h.c.} , (Bl) 
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where a{k) and a^{k) are the annihilation and creation operators of the scalar field, and p{k) is polynomially bounded 
in k. Introducing operators 

O^iio)^ J -^{gi.^ - ^k)p{k)a{k) ± g{u: + LO^)p{k)a\k)] , (B2) 

a calculation using the commutation relations satisfied by the annihilation and creation operators and the convolution 
theorem in the form 

-^g{u^ - J)g[u;') = g^u;) = J[u) (B3) 



(and the fact that g{u) — g{—u) as g is real) yields the expression 

Since the right-hand side of Eq. IIB4|) is (formally) a manifestly positive operator, we conclude that the expectation 
value of in any quantum state must satisfy the following bound: 

^^^^''^^r^ J (^^i^i'^+^opbwp (B5) 

for all (physically reasonable) states ip. It follows that the averaged energy density satisfies a state-independent lower 
bound. Further manipulation leads to the bound of Eq. I|l.lf) with 5—1. 
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